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Method of processing seismic data acquired by means of 
multicomponent sensors 

GENERAL FIELD 

5 The invention relates to techniques for processing 
. seismic data acquired by means of multicomponent 
sensors . 

This invention is in particular applicable to 
10 acquisition by means of cables disposed on the bottom 
of the sea (so-called "OBC" or "Ocean Bottom Cable" 
techniques) . 

Multicomponent geophones capable of working in any 
15 position whatsoever, in particular at the bottom of the 
sea, have recently been proposed. This "omnitilt" probe 
technology has allowed new simplified cables 
(mechanical joints are no longer necessary) and allows 
acquisitions with a better seismic bandwidth. 

20 

However, the acquisition step does not make it possible 
to provide the true orientation of the geophones of the 
cable, although this information is indispensable for 
making it possible to process the data. 

25 

The invention proposes a processing which is intended 
to be implemented on raw data and which allows 
reorientation and calibration (intended to convert the 
measurements of various geophones into a common phase 
30 and amplitude response) . 

STATE OF THE ART 

Techniques consisting in isolating from the signal the 
35 data which correspond to the first arrival at the 
sensor and in determining on the basis of these data a 
filter intended to be applied to the raw data so as to 
correct them and to thus obtain the components of the 
signal on the expected axes have already been proposed. 
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A proposal to this effect has been described in the 
article : 

"Horizontal vector infidelity correction by general 
5 linear transformation" - Joe Dellinger et al. - SEG - 
9-14 September 2001. 

However, this technique is not necessarily optimal 
since the coupling mechanism which intervenes at the 
10 geophone level is not the same for the waves which 
correspond to a first arrival at the sensor and for the 
waves reflected or converted by the seismic horizons. 

PRESENTATION OF THE INVENTION 

15 

The invention proposes another approach which employs 
the true data window for numerically reconstructing 
geophones oriented along the desired axes. 

20 Implicitly, this approach compensates for the errors 
which are not related to the geophones themselves, but 
which are due to the fact that the coupling between the 
geophone and the waves to be recorded is different 
depending on whether it is necessary to make a vertical 

25 vibration movement rather than a horizontal movement 
(on account of gravity) . 

In the case of a cable, the coupling is furthermore 
different depending on whether the vibration movement 
30 is in the direction of the cable or transverse. 

Moreover, since deeper windows are subject to a lower 
S/N (signal-to-noise) ratio, processing which 
implements trace stacks is moreover used. 

35 

The invention thus proposes, according to a first 
aspect, a method of processing seismic data acquired by 
means of a sensor having at least three geophone 
components, characterized in that estimators are 
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determined which are combinations of these components 
making it possible to isolate the various data 
depending on whether they correspond to propagation 
with reflection or with conversion and in that, to 
5 determine a sensor reconstruction, the operators to be 
applied to the various components of the sensor are 
determined in such a way as to minimize the deviation 
between reference data and data obtained by applying 
the estimators to the sensor reconstruction, the 
10 operators thus determined being applied to the data 
acquired. 

It is specified here that, in the remainder of the 
present text, the term geophone is understood to mean 
15 any velocity sensor and the term hydrophone any 
pressure sensor. 



Preferred, but non limiting aspects of the method 
according to the first aspect of the invention are the 
20 following: 

the sensor furthermore including a hydrophone, the 
reference data for reconstructing a vertical geophone 
are derived from the data acquired by the hydrophone; 

the reference data for reconstructing a vertical 
25 geophone without hydrophone or for reconstructing 
horizontal geophones are derived from the application 
of the estimators to one of the geophones of the 
sensor; 

the orientation in the horizontal plane of a 
30 geophone component is obtained by minimizing the 
estimator of the transverse reflection; 

the estimators are determined as a function of a 
model of isotropic propagation or including the 
azimuthal ani sot ropy . 

35 

According to another more general aspect, the invention 
proposes a method of processing seismic data acquired 
by means of a sensor having at least three geophone 
components, characterized in that estimators are 
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determined which are combinations of these components 
making it possible to isolate the various data 
depending on whether they correspond to propagation 
with reflection or with conversion. The estimators thus 
5 determined may find applications other than that 
forming the subject of the method according to the 
first aspect of the invention. 



DESCRIPTION OF THE FIGURES 

10 

Figure 1 is a diagrammatic representation giving 
the angular conventions used; 

- Figures 2 and 3 are flow charts giving the various 
steps of the processing respectively in one and the 
15 other of the two exemplary implementations described. 

DESCRIPTION OF ONE OR MORE MODES OF IMPLEMENTATION OF 
THE INVENTION 

20 First exemplary implementation : case of an isotropic 
propagation model 



Under the assumption of a locally ID (one dimensional) 
geology in proximity to the receivers, and assuming 
25 isotropic propagation of the earth, a given geophone, 
with an orientation $ measures: 



ntk = Rpp cos (\|/) 8pp k + (Rps cos (9 k - <{>) + Rtrsv sin(0 k - 
<|>) )sin(\|/)5ps,c 

30 

With : 

k : index for the shotpoint (from 1 to N) 
0 k : azimuth of the shotpoint with respect 

to the abscissa axis 
35 -Rpp. : reflectivity PP 

5 PP : dynamic correction PP ("normal 

moveout" or NMO) 
Rps ; isotropic radial reflectivity PS 
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Rtrsv : isotropic 

reflectivity PS 
5 Ps : dynamic correction 
moveout" or NMO) 



PS 



transverse 



("normal 



10 



This- model allows the evaluation of the reflectivity 
parameters from the set of traces tr k through simple 
processing of least squares comparison (ignoring \|/ to 
begin with) in the Fourier domain, leading to the 
following equations : 



N 



wc(<p) 
ws(<p) 



wc(<p) ws(<p) 
Sc2(<p) Scs(<p) 
Scs((p) Ss2(<p) 





f Rpp^ 




f Svx > 




Rps 




Shcx(<p) 




K Rtrsv j 




^Shsx((p) j 



Scalar quantities : 

15 

Sc(<|>) = 2 k cos (e k -(j)) Ss(«J>) = E k sin(9 k -<|>) 
Sc2(<|)) = £ k cos 2 (9 k -<|>) Ss2 (<|>) = Z k sin 2 (e k -<j)) 
Scs(<|)) = 2 k cos (e k -<|)) sin(6 k -<|>) 
N = Sc2 + Ss2 

20 

Wavelet quantities : 

wc(<|>) = Z k cos (9 k -(|>) 8ps k 8pp k _1 

25 ws (<J>) = 2 k sin (0 k -(|))5ps k 8pp k _1 

Trace stack quantities for geophone x : 

Svx = 2 k x k 8pp k _1 
30 Shcx(<j>) = Z k cos (9 k -<|>) x k 6ps k -1 
Shsx(<|>) = E k sin (9 k -<|)) x k 5ps k _1 

The solution of this linear system gives: 

Aiso Rpp cos (y) = W Svx + (Scs ws-Ss2 wc)Shcx + (Scs 
35 wc-Sc2 ws)Shsx 

Aiso Rps sin(\j/) = Kis cos (i|>) - Kic sin(<j>) 
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With: 

W=Sc2 Ss2~Scs 2 
5 Aiso = WN + wc(Scsws - Ss2 wc) + ws(Scswc Sclws) 

Kic = (Scs ws - Ss2 wc)s vx + (n Ss2 -wsws )shcx + (-N Scs + wswc )Shsx 
Kis = (Scs wc - Sc2 wi)Svx: + (-N Scs + wc w5 )shcx + (n Sc2 - wc w^jshsx 

This modelling allows evaluations taking account of the 
following properties: 
10 a. Rpp does not depend on (]), 

b. | Rps 1 2 + 1 Rtrsv | 2 does not of course depend on (|> 
either, 

c. Aiso is in practice rapidly steady over time and 
can be ignored for the calibration/orientation 

15 procedure, since it is common to all the geophones of 
one and the same receiver. 

Evaluations of dense shots 

20 Most of the OBC acquisitions are gleaned using a dense 
and regular grid of sources, which allows considerable 
simplification : 

Sc = Ss = 0, wc = ws=0 (symmetry of the sources with 
respect to the receivers) 
25 Scs = 0 Sc2 = Ss2 = N/2 (isotropic source distribution) 

Next, the exact solution can be obtained through the 
approximation : 
N Rpp cos (\|/) = Svx 
30 N Rps sin(v|/) = 2 Shcx(<|>) 

N Rtrsv sin(\|/) = 2 Shsx(<|>) 

This approximation leads to very simple calculations, 
not involving any wavelets, and can be applied 
35 immediately. 



Orientation of the geophones 
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Since Rtrsv does not exist physically, the minimization 
of the energy of Rtrsv leads to a trigonometric 
equation which gives the true orientation <j> ge o( + k *) : 

5 tan (2 <|> g eo) = 2 (£Ki Cl Kis,)/(£Kic : 2 -£Kis, 2 ) 

( (Emax-Emin) / (Emax+Emin) ) 1/2 gives a check on the 
quality of the reorientation. 

10 Moreover, if one wishes to find the orientation 
according to the first arrivals, it is possible to 
correct the said first arrivals so as to set them to 
one and the same arrival time, then to simplify kic and 
kis by replacing the wavelets wc and ws by the scalars 

15 sc and ss, by considering that the waves recorded 
horizontally are in fact the projection of the radial 
wave P, present on all the geophones since it is 
oblique . 

20 Geophone vertical composite calibration : 

With the geophones gi, g2, g3/ we construct a vertical 
composite geophone v, v = opx*gi + op 2 *g2 + op3*g3 (or 
comprising additional similar terms in the case where 
25 extra geophones are present in the receiver) where opi, 
op2/ op 3 are the filters of finite length and op u *g u 
represents the convolution of geophone g u with filter 
op u . 

30 such that: 

El = |XH - XV] 2 - IKic(v) ) 2 + |Kis(v) | 2 

The energy of the difference between XH (hydrophone 
35 after application of the geophone phantom, or 
cross-ghost hydrophone) and XV ,(the vertical composite 
geophone after application of the hydrophone phantom or 



WO 2005/040857 PCT/IB2004/003703 

- 8 - 

cross-ghost geophone) , (see for example in this regard 
the Applicant's Patent Application FR 2 743 896). 

E2 = |Rps(v)| 2 + |Rtrsv(v)| 2 
5 horizontal energy of the vertical composite, 

E = XE1 + (1-A,)E2 is a quadratic form in the 
coefficients of the filters and can be reduced to the 
minimum, thus giving a linear system to be solved. (X 
10 is a matching parameter, 0<=A,<=1, which favours either 
a greater adjustment to the reference hydrophone or a 
greater minimization of the shear energy) . 



In the case of terrestrial data, that is to say if 
15 there is no hydrophone available, it is possible to 
choose one of the geophones as reference and to replace 
the hydrophone by Rpp(g r ef)- 

Calibration in a horizontal arbitrary direction 

20 

With g = opi*gi + op 2 *g2 + op 3 *g 3 and cp g an arbitrary 

direction, 

We define: 

El - |Rps(g, cp g ) - Rps(ref, cp ref )| 2 , 
25 as being the energy of the difference between the 
evaluation of Rps of the arbitrary composite geophone 
and the evaluation of Rps of a reference geophone (in 
general the geophone oriented in the direction of the 
cable) . 

30 

E2 = |Rpp(g)| 2 + |Rtrsv(g, cp g )| 2 (the nonradial energy) 

E = AJEl + (1 - X)E2 allows the derivation of a 
composite horizontal geophone in the desired direction, 
35 having the- .same frequency response as the reference 
geophone, and with a minimum PP contamination. 
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Considering the cases (p re f = 0 and cp re f - tc/2, it is 
possible by simple trigonometric combination to 
generate the radial and transverse projections. 



Second exemplary implementation 
azimuthal propagation modelling 



case of anisotropic 



Let a be the direction of the natural fast propagation 
axis. The modelling of the measurement of the geophone 
10 becomes (using one or other of Rpsl and Rps2 the two 
images along the natural directions, or Rps and 8Rps 
defined by Rpsl = Rps + 8Rps, Rps2 = Rps - SRps) : 

mfc = Rpp cos(\|;)8pp k + (Rpsl cos (cp-a) cos (9 k -a) + Rps2 
15 sin (cp-a) sin (0 k -a) ) sin (v|/) 8ps k 

nik = Rpp cos (\j/) 8pp k + (Rps cos (0k~cp) + 8Rps cos (0 k +<p- 
2a) ) sin (\|/) Sps k 



20 



giving the normal equations 



M = 



N 



wc 



wc 



GO 



M. 



(<p) wc(2a-q>) 

Sc2(q>) Sc2(a)-Nsin 2 (a-cp) 
ws(2a-<p) Sc2(a)-Nsin 2 (a-<[>) Sc2(2a-cp) 

\ 



^Rpp^ 


f 


Rps 




^5Rps > 


V 



Svx 
Shcx(cp) 
Shcx(2a-q>) 



The solution of this linear system gives: 

Aiso Rpp cos (\|/) = unchanged 
25 Aaniso Rps sin(\j/) = (Kac cos (2a-(p) + Kas sin (2a- 
<p) ) sin (2 (a-cp) ) 

Aaniso Rtrsv sin(\|/) « (-Kac cos (cp) -Kas sin (<p) ) sin (2 (a- 
q>) ) 



30 With: 
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Aaniso « sin 2 (2 (oc-cp) ) Aiso 

Kac = (Sc2 ws-Scswc ) Svx + (N Scs-wc ws ) Shcx - (N Sc2- 
wcwc ) Shsx 

5 Kas = (-Ss2^+Scsvw? ) Svx + (N Ss2-ws ws ) Shcx - (N Scs- 
ws wc ) Shsx 

Vertical calibration of composite geophone : 

10 The isotropic process remains applicable with the 
change 

E2 = |Kac(g)| 2 + |Kas(g)| 2 

Horizontal arbitrary calibration or rows/columns of 
15 composite geophones 

The observation of SRps over the data field makes it 
possible to diagnose the presence (or otherwise) of 
significant azimuthal anisotropy, (the quantity 
20 sin 2 (2 (a-cp) ) SRps does not require a knowledge of a for 
its calculation) . 

The isotropic process remains applicable with the 
changes 

25 E2 = |5Rps(v)| 2 and E = X (E1+E2) + (1-X) E3 . 

When a is not generally known, a scan over a range of 
7i/2 is implemented, using the value of a which 
minimizes Emini/EO . 



